Obtaining Electronic Properties of Molecules through Combining Density Functional Tight Binding with Machine Learning.

We have introduced a machine learning workflow that allows for optimizing electronic properties in the density functional tight binding method (DFTB). The workflow allows for the optimization of electronic properties by generating two-center integrals, either by training basis function parameters directly or by training a spline model for the diatomic integrals, which are then used to build the Hamiltonian and the overlap matrices. Using our workflow, we have managed to obtain improved electronic properties, such as charge distributions, dipole moments, and approximated polarizabilities. While both machine learning approaches enabled us to improve on the electronic properties of the molecules as compared with existing DFTB parametrizations, only by training on the basis function parameters we were able to obtain consistent Hamiltonians and overlap matrices in the physically reasonable ranges or to improve on multiple electronic properties simultaneously.

Fano Resonance and Incoherent Interlayer Excitons in Molecular van der Waals Heterostructures.

Complex van der Waals heterostructures from layered molecular stacks are promising optoelectronic materials offering the means to efficient, modular charge separation and collection layers. The effect of stacking in the electrodynamics of such hybrid organic-inorganic two-dimensional materials remains largely unexplored, whereby molecular scale engineering could lead to advanced optical phenomena. For instance, tunable Fano engineering could make possible on-demand transparent conducting layers or photoactive elements, and passive cooling. We employ an adapted Gersten-Nitzan model and real time time-dependent density functional tight-binding to study the optoelectronics of self-assembled monolayers on graphene nanoribbons. We find Fano resonances that cause electromagnetic induced opacity and transparency and reveal an additional incoherent process leading to interlayer exciton formation with a characteristic charge transfer rate. These results showcase hybrid van der Waals heterostructures as paradigmatic 2D optoelectronic stacks, featuring tunable Fano optics and unconventional charge transfer channels.

Revealing generation, migration, and dissociation of electron-hole pairs and current emergence in an organic photovoltaic cell.

Using an innovative quantum mechanical method for an open quantum system, we observe in real time and space the generation, migration, and dissociation of electron-hole pairs, transport of electrons and holes, and current emergence in an organic photovoltaic cell. Ehrenfest dynamics is used to study photoexcitation of thiophene:fullerene stacks coupled with a time-dependent density functional tight-binding method. Our results display the generation of an electron-hole pair in the donor and its subsequent migration to the donor-acceptor interface. At the interface, electrons transfer from the lowest unoccupied molecular orbitals (LUMOs) of thiophenes to the second LUMOs of fullerene. Further migration of electrons and holes leads to the emergence of current. These findings support previous experimental evidence of coherent couplings between electronic and vibrational degrees of freedom and are expected to stimulate further work toward exploring the interplay between electron-hole pair (exciton) binding and vibronic coupling for charge separation and transport.

Coherent Real-Space Charge Transport Across a Donor-Acceptor Interface Mediated by Vibronic Couplings.

There is growing experimental and theoretical evidence that vibronic couplings, couplings between electronic and nuclear degrees of freedom, play a fundamental role in ultrafast excited-state dynamics in organic donor-acceptor hybrids. Whereas vibronic coupling has been shown to support charge separation at donor-acceptor interfaces, so far, little is known about its role in the real-space transport of charges in such systems. Here we theoretically study charge transport in thiophene:fullerene stacks using time-dependent density functional tight-binding theory combined with Ehrenfest molecular dynamics for open systems. Our results reveal coherent oscillations of the charge density between neighboring donor sites, persisting for ∼200 fs and promoting charge transport within the polymer stacks. At the donor-acceptor interface, vibronic wave packets are launched, propagating coherently over distances of more than 3 nm into the acceptor region. This supports previous experimental observations of long-range ballistic charge-carrier motion in organic photovoltaic systems and highlights the importance of vibronic coupling engineering as a concept for tailoring the functionality of hybrid organic devices.

Discriminative modulation of the highest occupied molecular orbital energies of graphene and carbon nanotubes induced by charging.

The highest occupied molecular orbital (HOMO) energies of carbon nanotubes (CNTs) and graphene are crucial in fundamental and applied research of carbon nanomaterials, and so their modulation is desired. Our first-principles calculations reveal that the HOMO energies of CNTs and graphene can both be raised by negatively charging, and that the rate of increase of the HOMO energy of a CNT is much greater and faster than that of graphene with the same number of C atoms. This discriminative modulation holds true regardless of the number of C atoms and the CNT type, and so is universal. This work provides a new opportunity to develop all-carbon devices with CNTs and graphene as different functional elements.

Time-dependent density functional theory quantum transport simulation in non-orthogonal basis.

Basing on the earlier works on the hierarchical equations of motion for quantum transport, we present in this paper a first principles scheme for time-dependent quantum transport by combining time-dependent density functional theory (TDDFT) and Keldysh's non-equilibrium Green's function formalism. This scheme is beyond the wide band limit approximation and is directly applicable to the case of non-orthogonal basis without the need of basis transformation. The overlap between the basis in the lead and the device region is treated properly by including it in the self-energy and it can be shown that this approach is equivalent to a lead-device orthogonalization. This scheme has been implemented at both TDDFT and density functional tight-binding level. Simulation results are presented to demonstrate our method and comparison with wide band limit approximation is made. Finally, the sparsity of the matrices and computational complexity of this method are analyzed.

Dissipative time-dependent quantum transport theory.

A dissipative time-dependent quantum transport theory is developed to treat the transient current through molecular or nanoscopic devices in presence of electron-phonon interaction. The dissipation via phonon is taken into account by introducing a self-energy for the electron-phonon coupling in addition to the self-energy caused by the electrodes. Based on this, a numerical method is proposed. For practical implementation, the lowest order expansion is employed for the weak electron-phonon coupling case and the wide-band limit approximation is adopted for device and electrodes coupling. The corresponding hierarchical equation of motion is derived, which leads to an efficient and accurate time-dependent treatment of inelastic effect on transport for the weak electron-phonon interaction. The resulting method is applied to a one-level model system and a gold wire described by tight-binding model to demonstrate its validity and the importance of electron-phonon interaction for the quantum transport. As it is based on the effective single-electron model, the method can be readily extended to time-dependent density functional theory.

Time-dependent density functional theory based Ehrenfest dynamics.

Time-dependent density functional theory based Ehrenfest dynamics with atom-centered basis functions is developed in present work. The equation of motion for electrons is formulated in terms of first-order reduced density matrix and an additional term arises due to the time-dependence of basis functions through their dependence on nuclear coordinates. This time-dependence of basis functions together with the imaginary part of density matrix leads to an additional term for nuclear force. The effects of the two additional terms are examined by studying the dynamics of H(2) and C(2)H(4), and it is concluded that the inclusion of these two terms is essential for correct electronic and nuclear dynamics.

Communication: Linear-expansion shooting techniques for accelerating self-consistent field convergence.

Based on the corrected Hohenberg-Kohn-Sham total energy density functional [Y. A. Zhang and Y. A. Wang, J. Chem. Phys. 130, 144116 (2009)], we have developed two linear-expansion shooting techniques (LIST)- direct LIST (LISTd) and indirect LIST (LISTi), to accelerate the convergence of self-consistent field (SCF) calculations. Case studies show that overall LISTi is the most robust and efficient algorithm for accelerating SCF convergence, whereas LISTd is advantageous in the early stage of an SCF process. More importantly, LISTi outperforms Pulay's direct inversion in the iterative subspace (DIIS) [P. Pulay, J. Comput. Chem. 3, 556 (1982)] and its two recent improvements, energy-DIIS [K. N. Kudin, G. E. Scuseria, and E. Cancès, J. Chem. Phys. 116, 8255 (2002)] and augmented Roothaan-Hall energy-DIIS [X. Hu and W. Yang, J. Chem. Phys. 132, 054109 (2010)].

Time-dependent density-functional theory/localized density matrix method for dynamic hyperpolarizability.

Time-dependent density-functional theory/localized density matrix method (TDDFT/LDM) was developed to calculate the excited state energy, absorption spectrum and dynamic polarizability. In the present work we generalize it to calculate the dynamic hyperpolarizabilities in both time and frequency domains. We show that in the frequency domain the 2n+1 rule can be derived readily and the dynamic hyperpolarizabilities are thus calculated efficiently. Although the time-domain TDDFT/LDM is time consuming, its implementation is straightforward because the evaluation of the derivatives of exchange-correlation potential with respect to electron density is avoided. Moreover, the time-domain method can be used to simulate higher order response which is very difficult to be calculated with the frequency-domain method.

Density matrix based time-dependent density functional theory and the solution of its linear response in real time domain.

A density matrix based time-dependent density functional theory is extended in the present work. Chebyshev expansion is introduced to propagate the linear response of the reduced single-electron density matrix upon the application of a time-domain delta-type external potential. The Chebyshev expansion method is more efficient and accurate than the previous fourth-order Runge-Kutta method and removes a numerical divergence problem. The discrete Fourier transformation and filter diagonalization of the first-order dipole moment are implemented to determine the excited state energies. It is found that the filter diagonalization leads to highly accurate values for the excited state energies. Finally, the density matrix based time-dependent density functional is generalized to calculate the energies of singlet-triplet excitations.